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FOREWORD 

This document presents the results of work performed 
by personnel of the Advanced Technology Systems Section of 
Lockheed's Huntsville Research & Engineering Center. This 
document was prepared for the Aerodynamic Systems Analysis 
Section of the NASA- Johnson Space Center, Houston, Texas un- 
der Contract NAS9-14517, Barney B. Roberts, technical moni- 
tor. 
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SUMMARY 

The governing equations for the one- dimensional flow 
of a gas-particle system are discussed* Gas-particle effects 
are coupled via the system momentum and energy equations 
with the gas assumed to be chemically frozen or in chemical 
equilibrium. A computer code for calculating the one-dimen- 
sional flow of a gas-particle system is discussed and a user's 
input guide presented. 

The computer code provides for the expansion of the gas- 
particle system from a specified starting velocity and nozzle in- 
let geometry. Though general in nature* the final output of the 
code is a startline for initiating the solution of a supersonic gas- 
particle system in rocket nozzles* The startline includes gasdy- 
namic data defining gaseous startline points from the nozzle cen- 
terline to the nozzle wall and particle properties at points along 
the gaseous startline. 
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NOMENCLATURE 


Symbol 

A 



m 

P 

Pr 

r 

R 

Re 

S 

T 


u 

V 

w 

X 


Description 

local area 
drag parameter 
heat transfer coefficient 
slip flow velocity factor 
drag coefficient 

specific heat at constant pressure 

total enthalpy 

static enthalpy 

mass 

pressure 

Prandtl number 

radius 

gas constant 

Reynold* s number 

entropy 

temperature 

velocity 

velocity of sound 
mass flow rate 

ratio of total particle mass flow 
rate to total gas mass flow rate 
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Symbol 


Description 


x 

V 


axial coordinate 

ratio of j ** 1 particle mass flow 
rate to total particle mass flow 
rate 


Greek 


y 


p 

Subscripts 

P 


m 


viscosity 

density 


particle 

mixture 


g 

r 

x 


Superscripts 

j 


gas 

radial component 
axial component 


.th 

J 


particle 


np 


number of particles 
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Section 1 
INTRODUCTION 

Solid propellant formulations are used in the Space Shuttle booster and 
stage separation motors. Interaction of the solid propellant exhaust plumes 
with the external surroundings can result in a hostile environment which must 
be considered in the sign process. The solid propellant boost motor exhaust 
significantly affects the vehicle base environment. During the separation of the 
orbiter and booster motors the exhaust plumes from the separation motors im- 
pinge on the orbiter. Particle impingement results in erosion of the surfaces 
as well as high heating rates and forces acting on the surfaces. Consequently, 
gas-particle flows in nozzle exhaust plumes have received considerable at- 
tention. 

A recent formulation (Ref. 1) has extended the supersonic two-dimen- 
sion, . axisymmetric calculation (Ref. 2) to include the treatment of chemically 
reacting flows. A computer code was subsequently developed which treated 
chemical equilibrium (Ref. 3) or chemical kinetics (Ref. 4) in the gas phase. 
Calculations are initiated along some data surface which must be everywhere 
supersonic. Gas data (velocity, flow deflection angle, pressure, etc.) and par- 
ticle data (velocity, density, streamline location, etc.) are required to initiate 
the solution. 

Originally, a computer code developed by Kliegel (Ref. 2) was used to 
calculate a startline to initiate a supersonic method-of-characteristics solu- 
tion. The Kliegel code utilizes a Sauer (Ref. 4) approximation to the transonic 
solution. One- dimensional particle lags (velocity and temperature) which are 
assumed constant from the chamber are determined. Particle limiting stream- 
line trajectories are then traced from some starting location upstream of tne 
startline to determine their intersections with the startline. The location of 
the startline is specified by input data. These data are then utilized in the 
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transonic approximation to obtain the gas-particle startline conditions. This 
calculation procedure was incorporated into the gas-particle reacting gas code 
(Ref. 5) to provide startline information. The Kliegel calculation assumes con- 
stant thermodynamic properties in the gas and particle phases. Combining the 
transonic data with chemical equilibrium thermodynamics (Ref. 3) results in a 
mismatch of thermodynamic properties and at least one state variable. To min- 
imize the mismatch, the gas velocity, pressure and temperature from the tran- 
sonic approximation and chemical equilibrium thermodynamic data are used to 
initiate the supersonic flow solution. The supersonic computer code calculates 
the change in entropy and total enthalpy levels. 

Solutions to the transonic approximation become difficult for throat ra- 
dius of curvature ratios less then 1.5. This difficulty, coupled with the constant 
thermodynamic property assumption, prompted the development of a reacting 
gas code (Appendix B) to provide startline data for the RAMP code (Ref. 5). 

The solution is one-dimensional and is intended for interim use until a fully 
coupled transonic solution for two-dimensional chemically reacting flows is 
developed. 

The one-dimensional gas-particle solution is fully coupled through the 
mixture momentum and energy equations. Particle distributions are repre- 
sented as a series of discrete sizes and/'"'** chemical species. The gas is as- 
sumed to be either in chemical equilibrium, or frozen or to have constant ther- 
modynamic and transport properties. The solution is initiated at a starting con- 
dition specified geometrically and gasdynamically. Analytical functions are used 
to describe the geometry from which the variation in area ratio ir» obtained. The 
work performed by the gas in accelerating the particles and the heat given up by 
the particles are a function of the relative distance over which the expansion oc- 
curs. The gas-particle flow is defined at specified area ratios and axial sta- 
tions relative to the initial data plane. At each axial station the mixture pro- 
perties are solved iteratively. The solution is continued until a gas Mach num- 
ber specified by the user is reached. Particle trajectories are then traced 
through the region of interest by numerically integrating the particle equations 
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of motion. The complete solution consists of the one-dimensional gas-particle 
properties and particle limiting streamline locations at a specified data line. 

The following sections describe the flow assumptions, present a develop- 
ment of the governing equations and outlines the user’s guide for a computer 
code (ODPART) to calculate one-dimensional, chemically reacting gas-particle 
flows. 
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Section 2 

TECHNICAL DISCUSSION 

The one-dimensional gas-particle solution is fully coupled via the mo- 
mentum and energy equations. Particle distributions are represented as a 
series of discrete sizes and/or chemical species. The gas is assumed to be 
in chemical equilibrium, frozen or to have constant thermodynamic and trans- 
port properties with no mass exchanged between the phases. The solution is 
initiated at a specified condition. Analytical functions are used to describe 
the geometry from which the variation in area ratio is obtained. The work 
performed by the gas in accelerating the particles and the heat given up by 
the particles are a function of the relative distance over which the expansion 
occurs. Thus the gas-particle flow is defined at specified area ratios and ax- 
ial stations relative to the initial data plane. At each axial station the mixture 
properties are solved iteratively. The solution is continued until a gas Mach 
number specified by the user is reached. Particle trajectories are then tra- 
ced through the region of interest by numerically integrating the particle equa 
tion of motion. The complete solution consists of the one-dimensional gas- 
particle properties and particle limiting streamline locations at a specified 
data surface. 

The following subsections describe the flow assumptions, present a de- 
velopment of the governing equatio and discuss a comparison of analytical 
calculations generated with start lines from the Kliegel transonic analysis 
and the ODPART compute 1 * code. 

2. 1 EQUATIONS DEFINING THE ONE -DIMENSIONAL FLOW OF A 
GAS-PARTICLE MIXTURE 

The flow of a gas -particle mixture is described by the equations for con 
servation of momentum, conservation of energy and conservation of mass. 
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Development of the equations defining the gas and particle phases and their 
coupling in the conservation equations is based on the following assumptions: 

1. The particles are spherical in shape. 

2. The particle internal temperature Is uniform. 

3. The gas and particles exchange thermal energy 
by convection. 

5. The gas obeys the perfect gas law and is either 
frozen or in chemical equilibrium. 

5. The forces acting on the control volume are the 
pressure of the gas and the drag of the particles. 

6. The gas is inviscid except for the drag it exerts 
on the particles. 

7. There are no particle interactions, i.e. , no col- 
lisions and no exchange of thermal energy. 

8. The volume occupied by the particles is negligible. 

9. There is no mass exchange between the phases. 

10. The particles are inert. 

11. A discrete number of particles, each of different 
size or chemical species, are chosen to repre- 
sent the actual continuous particle distribution. 


In a gas -particle flow, the coupling between the particles and gas is im- 
portant and must be properly treated in order to construct an accurate model. 
Particles in the flow gain motion from the drag forces exerted by the combus- 
tion gases and, transfer heat to the gas primarily by convection. Losses in a 
gas -particle flow result from the difference between particle and gas velocity, 
and from the difference between particle and gas temperature. Coupling be- 
tween the gas and particle phases in this analysis is accomplished through the 
momentum and energy equations. The one -dimensional momentum equation 
for the gas -particle mixture is 


pu 


du 

dx 



j j du 
p u 


) 


dP 

dx 


0 


(1) 
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The one-dimensional energy equation for the mixture is given by 


2 

np 


| h + 

+ £. p Va 

h j + u j 2 

i 2 

j=l 

2 


( 2 ) 


.ne remaining equation used to define the mixture properties is the conserva 
tion of mass equat ; on for a gas-particle flow 

np 

w = puA + ^ p^u^A (3) 

m frf 


Thermodynamic relations used in the solution of gas phase properties include 


dP 

TdS = dh - -Sf 
P 


(4) 


and pressure as a function of density and entropy 


dP = v 2 dp + q~~T~r dS 


( 5 ) 


For accurate modeling of gas-particle flows, particle properties must 
be 'mown as well as gas phase and mixture properties. To obtain a particle 
drag coefficient, C^, a -able look-up is done to obtain the ratio 

Stoke 8 

as a function of local Reynolds number. The resulting ratio is then used to 
calculate the particle drag coefficient using the following relationship to cor- 
rect for slip flow 


c j = / C D V Re j [(l+9.45r j )(H2.52g j ) + 3.03r j2 ] 

° V C D Stoke ^/ [(l+9.45rj)(l+3.78c j )+ (. 91 1)(4. + i 1.34T j )]r j 2 


( 6 ) 
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where the parameter, 7, is calculated as follows 



(u - u p-* ) 
Re j <5T 


( 7 ) 


In Eq. (6), the Stokes drag coefficient is used as a reference drag coefficient 
calculated by 


C J 

Stokes 



( 8 ) 


The particle drag coefficient enters the calculation of particle properties 
through the drag parameter, a, which is defined as 

a? = 4* 5 i/j C J m ^ r ^ (9) 

D p p 

Change in partic 1 ' ocity with respect to change in axial coordinate is de- 
termined by 


du 


j = 


(u - u^)dx 

U j 


(10) 


Heat energy exchange between the particle and gas phases is established by 
the particle l.ec.t transfer coefficient, , defined by 


J - —R 


c J (l + 0.2295 Re 


j(0. 55 Pr J°‘ 33 ) 


Pr J CjJ 


( 11 ) 


This relationship assumes that heat energy exchange between the gas and par- 
ticles occurs only by convection. A table look-up is performed for Prandtl 

number Pr, specific heat at constant pressure, c , and viscosity, V, as a func- 

P 

turn of gas velocity. Change in particle enthalpy with respect to change in ax- 
ial coordinate is calculated by 
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(U) 


dh j , 2 aid grlte 

3 J 


u 


where the thermal lag, AT** is defined by 


AT J = T J 
P 


(13) 


As can be observed, the particle size, r \ enters the calculation of both du^ 

P 

and dh** in Eqs. (10) and (12), respectively. The particle size then appears in 
the mixture momentum and energy equations (Eqs. (1) and (2)) through the par- 
ticle velocity and enthalpy terms. Particle size thus has a significant effect 
on f he coupling between the gas and particle phases and the resulting perfor- 
mance of a particular nozzle-propellant combination. 


From the basic equations defining the gas -particle mixture, the gas 
phase properties and the particle properties, expressions were derived to cal- 
culate gas and particle properties at successive axial locations in a rocket noz- 
z*' The expressions were derived in differential form in order to be compa- 
tible with the forward marching integration technique used in the computer so- 
lution of the equations. 


In the following derivations 'to simplify the equations), it is assumed 
that there is only one particle specie present thus eliminating the need for the 


* 

summation symbolism, / . , on all particle properties. Complete symbolism 

1 = 1 


tor multiple particTe species is used for the basic governing Eqs. (1) through (12), 
and for the difference forn >f the derived expressions presented in Appendix A. 


Beginning with the gas phase thermodynamic relationship defined by Eq, (5), 
an expression was derived for calculating gas velocity. 
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Pressure may be written as a function of density and entropy as follows 


dP = v 2 dp + - dS (14) 

°P " 

Solving for dS, Eq. (14) becomes 

c - R y 

dS = p p ■ (dP - v dp) 

Multiplying through by temperature, T, and using the perfect gas law, P =pRT, 
the relation becomes 


TdS=-^r (c -R) (dP- v 2 dp) 

pK p 

Substituting for TdS using the basic thermodynamic relation TdS = dh - dP/ p, 
the equation is 

dh ~ = tW k - R > (dp • yZ > 

p pR p 

Dividing through by dx and rearranging terms the expression becomes 


dh _ dP = _L {c 
dx p dx pR * p 


ni / dP 2 dp , 

R > <sr • v $> 


(15) 


The energy equation for a gas-particle mixture is 


2 * j . j 2 

H = (h + ^) + V (h j +4-) 
£ w * 


(16) 
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Differentiating and dividing through by dx, the energy equation becomes 


dh + t (fli* * o 

dx dx w dx dx 


( 17 ) 


The momentum equation for a gas -particle mixture is 


du . i j du J . dp n 

/>»*;*?* ’ ar + & 0 


08) 


Rearranging terms, the momentum equation becomes 


du __ _1_ dP _ j du ^ 

u dx ~ p dx p dx 


(19) 


The ratio of particle to gas mass flow is written 


p^u - * A w^ 


pu A 


Rearranging terms and solving for the ratio o'" densities, the relation is 


j • j 

p J u w J 


( 20 ) 


Substituting Eq. (20) in Eq # (19), the momentum equation becomes 


u 


du _ _1 dP _ ^ j d\r^ 


dx 


p dx 


u 


1 '%v 

J w 


dx 


( 21 ) 
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Simplifying momentum Eq. (21) and substituting Eq. (21) in Eq. (17), the energy 
equation becomes 

. i - ±1„ *ti + «i(dbi + u j dui) =0 (22) 

dx p dx $, dx w 'dx dx / 

Collecting terms and simplifying, the equation is written 

4 - i ffl » j ^ - u £ii) = 0 (23) 

dx p dx ty \dx dx dx / 

Substituting Eq. (15) in Eq. (23), the energy equation becomes 

J_ (c . R) /^. v 2 d£v + f^/d^ + u j^ . u du^\ = Q 
pR ' p '\dx dx* vv 'dx dx dx f 


Collecting terms and multiplying through by R/(Cp - R) the relation is written 


J. dP . v^dp ( / R t w^ /dh^ 
p dx pdx \c - R / J, \ dx 


“ J #- • Sr) ■ ° 


Introduce the conservation of tnass relationship for the gas 


pu A = w 


and implicitly differentiate the equation 


uAdp + pAdu+pudA=0 


Dividing through by puAdx and rearranging terms, the expression becomes 


I d£ /I du , 1 dAl 
pdx - \u dx A dx / 
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Substituting Eq. (25) in Eq. (24), the energy equation is defined by 


1 dP 2 

- v 


p dx 




Solving for dP/pdx, the expression is 

2 j • 2 

Oft / 

dx ~il 

P 

Prom Eq. (21) the expression for dP/pdx is substitu ed in Eq. (26) 


idP v 2 du v 2 dA / R \ /w j \ /db j . jdu j du j \ .... 

p dx” u dx A dx '(c -R/( \dx + u dx ’ “dx ) * 26 * 


.. du u w J _j du J v 2 du v 2 dA 
U j " - U " = ** <■ " "T »' “ 


dx u j * 


OX 


u dx A dx 


/ R W&’k / dh ^ . j du ^ du 1 ^ 

(dT tuJ E- - u ar 


Collecting terms and rearranging, the relationship is 

du * 2 du v 2 u 2 dA wif / R \ / dh^ . j du ^ du^ \ du ^ 1 
dx "u 2 dx_ Au 2 ^ " w [( c p :R)( dx IT - u Hr)‘ u drJ 

Substitutions are made using the definition of Mach number, I' = u/v 


a 0 


du / 

. 2 
_ i . .a 

dA 

* j 
w J 

/ R \fdh j j 

du^ 

dJ) 

du ^ 

dx V 

M 2 ' AM 2 

dx 

• 

w 

l c p " R/\dx +u 

dx 

dx / 

dx 


= 0 


Substituting the basic relations for changes in particle velocity, d\. the re- 
lationship is written 


/ M 2 - 1 \ 

2 

u 

dA 


/ R \/ 2 

' M 2 f 

AM 2 

dx 

* 

w 

( c p - R/ \ 3 


- u a 


Jin* 


u 


j*) 


a a 


j (ujjdi 

.J 


(27) 
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Define a parameter, B 1 *, 


B- 


= ( c ) ( \ ( T '*' T ) * ( u " u ^) + u ( u * u ^)) + u(u-uj) 


Substituting B^ is Eq. (27) and multiplying through by M , the equation becomes 


du ,„,2 ,, M 2 u 2 dA , *» i*,2 _ rt 

^ ‘am 2 ” n 


Dividing through by u and multiplying through by dx, the relation is defined by 

, .2 . , u , A . B J M^dx 

du(M - 1) - dA + = 0 


Solving for du. 


da = (a dA - - 4 — (28) 

\ A U / - l 


an equation is developed from which gas velocity is calculated for the gas- 
particle mixture. 

Starting with the momentum equation an expression was derived for cal- 
culating gas pressure for the gas-particle mixture 

pudu + p-u J du^ + dP = 0 

Multiplying through by A, the momentum equation is 

puAdu + Ap^u^du^ + AdP = 0 (29) 
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Multiply differentials of the conservation of mass equations for the gas and 
particle phases by u and u^ respectively 

ud(puA) = 0 (30) 

u j d(p j u j A) = 0 (31) 

Equations (29), (30) and (31) are added to obtain 

puAdu + ud(puA) + p^u** Adu^ + u^d(p^u-^A) + AdP « 0 

Expanding the derivatives, d(puA) and d(p J u^A) the equation becomes 

puAdu + u Adp + puAdu + pu dA + p^u^Adu^ + u 1 * Adp^ + p^u^Adu^ + p^u^ dA 

+ AdP = 0 

Rearranging and collecting terms the expression is 

2puAdu + u^Adp + pu^dA + 2p^u^Adu^ + u^Adp + u^p^dA + AdP = 0 (32) 

The differential of the product PA is written 

d(PA) = PdA + AdP 

Solving this expression for AdP* the relation is written 

AdP = d(PA) - PdA (33) 

14 
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2 

The differential of pu A is written 

d(pu 2 A) = u 2 Adp + 2puAdu + pu 2 dA (34) 

* 2 

and the differential of p^u^ A is written 

d(p^u^A) = u^ 2 Adp + 2p^u^Adu^ + p^u^dA (35) 

Substituting Eqs. (33), (34) and (35) in Eq. (32), Eq. (32) becomes 
d(pu 2 A) + d(p J u j 2 A) + d(PA) - PdA = 0 


From the conservation of mass relationship the expression p^ = puwVu^w is 
obtained and substituted in the previous equation 

d(pu 2 A) + d . A) . + d(PA) - PdA = Q (36) 

u J w 


np . np t . 

Using the relationship, w Vw = X and wV w* Eq. (36) becomes 

j=i 3=1 


np 


d(pu A) + X ^ v J (o*v ’A) + d(PA) - PdA = 0 


The equation used to calculate gas pressure is then written 


d(PA) = PdA - d(pu A) - ^ ^v J (puu J A) 

3 = 1 


(37) 
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A relationship to calculate gas density was derived as follows from the 
conservation of mass equation for the gas 

puA = w 

Differentiating implicitly and dividing through by puA, the conservation equa- 
tion becomes 


d£ + du + dA _ Q 
p u A 

Solving for the change in density, a relationship is obtained for calculating gas 
density at successive integration steps 

<*> *-'(*♦ t) l38 > 


Gas temperature is calculated frora the perfect gas law using the pre- 
viously calculated values of density and pressure 

T = ^ (39) 

Entropy is calculated using the basic thermodynamic relationship de- 
fined by Lq. (5) 


dP = v 2 dp + dS 

C " A 

P 

Solving for dS, an equation is obtained for calculating entropy at successive 
integration steps 


dS 


c - R 

-Bp— (dP - v^dp) 


( 40 ) 
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Gas static enthalpy is calculated using the thermodynamic relationship 
defined by Eq. (4) 


TdS = dh 



Solving for dh, an equation is obtained for calculating gas static enthalpy from 
previously calculated values of pressure, density, temperature and entropy 

dP 

dh = TdS + — (41) 

P 


The gas total enthalpy is calculated by adding the gas static enthalpy to 
the gas velocity contribution 


2 

H = h + (42) 

2.2 EQUATIONS DEFINING THE CALCULATION OF PARTICLE LIMITING 
STREAMLINES 

If the user of the ODPART computer code flags the appropriate option, 
particle limiting streamline calculations are performed after the one-dimen- 
sional gas-particle calculations. The equations governing the particle limiting 
streamline calculations were obtained from the derivations of Ref. 1. The fol- 
lowing equations define compatibility relations along particle limiting stream- 
lines. Change in particle radial location along a limiting streamline is calcu- 
lated as a function of x using the following relationship 

dr i = d*jHi£ (43) 

U^X 
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A change in particle axial velocity as a function of x is given by 


du ^ = dx^ a J 
x 


(U 


- U J ) 
X 


u 


( 44 ) 


The change in particle radial velocity along a limiting streamline as a function 
of x is 


du 


j 

r 


= dx 


J a J 


(u - U-J.) 
r r 


u * 


(45) 


Change in particle enthalpy is given by 


dh^dx’ aV llillL 

u J J 


(46) 


The difference form of equations (43) and (46) used in the one- dimensioned com- 
puter code are presented in Appendix A. The particle equations of motion de- 
fined in Eqs. (43) through (46) are integrated from the initial data line to the 
calculated startline. The equations are solved iteratively at each social loca- 
tion stored during the one- dimensional gas-particle mixture calculations. 

The particle limiting streamline concept implies that there is a stream- 
line above which particles larger than a certain size cannot flow. The concept 
of a nozzle flowing full implies that particles of all sizes are distributed from 
the nozzle centerline to the nozzle wall. Both concepts of particle distribu- 
tions in two-phase flows are incorporated in the one-dimensionsd computer 
code presented in Appendix B. The user of the code may select the particle 
distribution scheme appropriate to a given application. 
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2.3 RESULTS OF ANALYTICAL COMPARISONS 


To confirm the validity of the one -dimensional, chemically reacting gas- 
particle calculations, an analysis was conducted to compare results obtained 
from the RAMP supersonic code using startlines generated with the Kliegel and 
ODPART computer codes. Three problem cases were analyzed with RAMP su- 
personic nozzle calculations generated with Kliegel and ODPART startlines for 
each case. The first case comparison used full-scale Space Shuttle solid rocket 
motor (SRM) geometry, ideal gas (constant property) thermodynamics, a cham- 
ber pressure of 500 psia and an aluminum propellant loading of 16%. The se- 
cond case comparison used full scale SRM geometry, real gas (chemically re- 
acting) thermodynamics, a chamber pressure of 500 psia and an aluminum pro- 
pellant loading of 16%. The final case comparison used a subscale nozzle, real 
gas thermodynamics, a chamber pressure of 1000 psia and an aluminum propel- 
lant loading of 2%. In each case, aluminum oxide was the particle species repre- 
sented by a six particle size distribution. A detailed report of the analysis is 
presented in Ref. 6. 

For ease of reference, the results of one case comparison are presented 
in Figs. 1 through 5. The results presented were calculated for full scale SRM 
nozzle geometry, real gas thermodynamics, a chamber pressure of 500 psia 
and a 16% aluminum propellant loading. The aluminum oxide particle species 
was represented by a six-particle size distribution with sizes ranging from a 
radius of 3.15 microns to 9.70 microns. 

In Fig. 1 the non-dimensional nozzle wall static pressure distributions 
for the RAMP nozzle calculations are compared. The RAMP calculation initi- 
ated with an ODPART startline demonstrates good agreement with the RAMP 
c alculation initiated with a Kliegel startline. The nozzle centerline static pres- 
sure distributions for the RAMP calculations are in excellent agreement in Fig. 

2. Location of particle limiting streamlines for the largest (9.70 micron radius) 
and smallest (3.15 micron radius) particles are compared in Fig. 3. Particle 
velocity and temperature distributions along the nozzle centerline are compared 
in Figs. 4 and 5, respectively. Agreement of particle properties between RAMP 
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calculations initiated with Kliegel and ODPART startlines is acceptable. Based 
on these comparisons it was concluded that the one -dimensional gas-particle 
calculations produce acceptable agreement with current analytical techniques. 
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Section 3 

CONCLUDING REMARKS 

An analytical technique has been developed for calculating the one-dimen- 
sional :low of a gas-particle system, '"he technique employs a fully coupled so- 
le with particle effects entering the calculations via the system momentum 
an rg^ equations. The one-dimensional equations are not re. ricted by noz- 
zle throat rad us of curvature rado as are current two-din nsional transonic 
calculations. The gas may be assumed to be chemically reacting, thus elimin- 
ating the problem of combining thermodynamic data from chemically frozen tran- 
sonic calculations with chemically reacting supersonic flow solutions. Particle 
lags (velocity and temperature) are allowed to vary, thus eliminating the con- 
stant lag assumption present in currently used transonic calculations. 

T1 e computer code (Appendix B) developed using this analytical approach 
expands a gas-particle system from a specified starting velocity and inlet geom- 
etry. Altkougb the code may be used to calculate the general one- dimensional 
flow of a gas-particle system, the final output of the code is a startline for in- 
itiating the solution of a supersonic gas-particle flow. The startline includes 
gasdynamic data defining gaseous startline points from the nozzle centerLue 
to the nozzle wall and particle properties at points along the gaseous startline. 

Results of supersonic nozzle calculations initiated with startlines from the 
ODPART code compare favorably with those initiated with startlines from the 
Kliegel transonic approximation. Based on these analytical comparisons (Ref. 6), 
it was concluded that the one-dimensional chemically reacting gas-particle sol- 
ution produces a valid and useful solution to two-phase flows in the inlet and 
throat regions of rocket nozzles. 
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DIFFERENCE EQUATIONS FOR CALCULATING THE FLOW OF 
A ONE-DIMENSIONAL GAS-PARTICLE MIXTURE 
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Appendix A 

The following equations are the difference form of particle Eqs. (10) and 
(12) presented in the main text of this report and of gas Eqs. (28), (37), (38), (40) 
and (41) derived in the text. The following equations are coded in the ODPART 
Computer code to calculate gas and particle properties at successive integra- 
tion steps. 


Gas Velocity : 


Au 


u . . . M 

— AA - dx 

A u 




1 

M 2 -l 


Gas Pressure: 


AP = - pu I Au + X 



V J a J 


(u-u j ) l 

u j J 


Gas Density: 


A p 




AA 
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Gas Static Enthalpy: 


Ah = - uAu + Ax X 
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Ga»? Entropy: 


AS ^ ^ (Ah 



Particle Velocity: 


Au 


j 



Ax 


Particle Enthalpy : 


Ah j = -|a j c j a,, 

3 u J 


Particle Limiting Trajectories: 


Ar^ = 


Ax^ 
r 
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Au J 
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Ax^** (u - u* 1 ) 

X X 
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Au J 
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Ax J a J (u r - u J r ) 


U‘ 


Ah'* = Ax* 


if 2 a J c*(T* - T) 
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Appendix B 

ODPART - A ONE DIMENSIONAL GAS-PARTICLE 
COMPUTER CODE 
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Appendix B 

B. 1 DESCRIPTION OF THE COMPUTER CODE 

ODPART is a computer code written in FORTRAN V language which per- 
forms one-dimensional calculations of gas-particle flows. The code was writ- 
ten for use on the Univac 1108 Executive VIII multiprocessor system but is eas- 
ily adapted for use on other systems. 

• Program Operation 

A schematic diagram of the ODPART solution is presented in Fig. A-l. 
The input data are read from cards and the gas and particle properties on the 
initial data line calculated. The gas and particles are assumed to be in dynam- 
ic and thermal equilibrium on the initial data line. A forward marching inte- 
gration scheme is employed for integrating the fully coupled equations defin- 
ing the gas-particle system. Solution of the equations is accomplished at suc- 
cessive axial stations at intervals equal to the input parameter, DX. At each 
new axial station, the first iteration is initialized with the properties of the pre- 
vious data line. The particle and gas properties are calculated and a check is 
made for convergence of the properties. If the solution has not converged, the 
program iterates the solution and recalculates the particle and gas properties 
based on the data from the prsvlov- iteration. If a solution at an axial station 
is not reached in 100 iterations., a*, a 1 or message is printed and » u e run is ter- 
minated. When the solution at a:< axial station converges, the converged values 
of the gas and particle properties are stored and output on a line printer. A 
c.ieck is made to determine whether or not the solution has progressed to the 
terminating Mach number specified by the programer. If the specified Mach 
number has not been reached, the solution is initialized for tne next axial sta- 
tion. When the specified Mach number is reached, the one-dimensional flow 
field calculations are terminated. By the appropriate choice of input flags, the 
user may instruct the program to calculate particle trajectories through the 
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flow field or distribute the particles along a gaseous start line from the nozzle 
centerline to the nozzle wall. If ICON(l) is set at 1, the program omits the par- 
ticle trajectory tracing calculations and the nozzle is assumed to be flowing full. 
A gaseous start line is constructed based on the value of ICON(12). If ICON (12) 
is 0. the program distributes gas points on the start line by a sine distribution. 

If ICON(12) ?s 1, the gas points are distributed at even radial increments. With 
ICON(l) as 1, particle start line data points are calculated for each gas point. 

The angle of the particle trajectory at each point is set equal to the gas flow 
angle. The gas and particle start line points are arranged in order and output 
in units compatible with the RAMP code input. If ICON(ll) is 0, the gas and 
particle start line points are output on a line printer and a card punch. If 
ICON(ll) is 1. the start line points are output on the line printer only. 

When ICON(l) is 0, the flowfield calculations are followed by particle tra- 
jectory tracing calculations. Particle trajectories are traced through the one- 
dimensional flowfield data previously calculated and stored. Particle trajectory 
calculations are performed from the initial axial station to the final axial sta- 
tion. The intersection of the particle trajectories with the line defining the fi- 
nal axial station provide the necessary data points to construct a start line of 
particle properties. As outlined previously a gaseous start line is constructed 
with gas data points either distributed at even radial increments or with a sine 
distribution. Each gaseous start line point consists of a radial co-ordinate, ax- 
ial co-ordinate, Mach number, flow angle, entropy, Mach angle and total enthalpy 
arranged in the necessary order and output with the proper units to be compati- 
ble with the input of the RAMP computer code. The intersections of the particle 
trajectories with the gaseous start line determine the location of particle limi- 
ting streamline points on the start line. If a particle limiting streamline inter- 
sects the start line between gas points, a new gas point is calculated at the in- 
tersection by interpolating between the two existing gas points. The new gas 
point is added to the gaseous start line. For the particle start line data points 
up to the limiting streamline intersections, particle properties are calculated 
at the axial and radial location of eacn gas point. Particle start line properties 
output for the equilibrium version of the RAMP code are axial and radial velocity 
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components, particle trajectory angle and particle density. Properties output 
for the finite rate version of the RAMP code are particle enthalpy and density, 
ratio of particle velocity to gas velocity and particle trajectory angle. As out- 
lined previously the gas and particle start line data points are output on a line 
printer or a line printer and caret punch depending on the value of ICON(ll). 

• Program Input Data 

Input data required for the calculations include program control flags and 
starting parameters, problem geometry, gas and mixture thermodynamics, and 
particle properties. Program control flags control the input to be read, the in- 
terval at which calculations are printed and stored, the axial distance between 
stations, and the type of start line which is output. The input starting para- 
meters determine the axial location of the initial data line relative to the noz- 
zle throat and the gas and particle velocities on the initial data line. 

The input thermodynamic data consist of thermodynamic and transport 
property data as a function of Mach number for the gas-particle mixture and for 
the gas phase. The gas-particle mixture thermodynamics includes Mach num- 
ber, molecular weight, ratio of specific heats, temperature, pressure, Prandtl 
number, viscosity and specific heat at constant pressure. Gas phase thermo- 
dynamics includes Mach number, gas constant, ratio of specific heats, Prandtl 
number, viscosity and specific heat at constant pressure. The user may chose 
to use either ideal gas thermodynamic relationships or real gas thermodynamic 
data. When ideal gas thermodynamics is s reified, a temperature exponent is 
input in place of the specific heat at constant pressure in the gas phase tnermo- 
dynamics. For ideal gas thermodynamics, only one entry is made in the gas- 
particle mixture and gas phase thermodynamic tables. Thermodynamic pro- 
perties in each table are input for a Macn number of 0.0. If read gas thermo- 
dynamics are used, one table of Mach numbers and corresponding properties 
is input for the gas -particle mixture thermodynamics and one or more tables 
are input for the gas phase thermodynamics. A maximum of 13 Mach number 
entries is permitted in each table. 
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The work required to accelerate the particles combined with the heat ex- 
changed between the particle and ^,as phase is a non isentropic process* In gas- 
particle nozzle expansions, the gas total enthalpy and entropy vary. To account 
for this variance analytically, tables of these variables may be generated by ex- 
panding the gaseous products from various reference conditions (total enthalpy 
and entropy). The ODPART code has the capability to interpolate for gas con- 
stant, ratio of specific heats and transport properties between themodynamic 
tables of data expanded from different chamber enthalpies and entropies. A 
maximum of ten tables with different total enthalpies may be input for the gas 
phase thermodynamics. For each total enthalpy, a maximum of two tables with 
different values of entropy are permitted. ODPART was written to accept the 
real gas thermodyanmic card output of the CEC computer code (Ref. 8). ODPART 
does not read a magnetic tape for thermodynamic data. 

The geometry defining a given problem establishes the boundary of the one- 
dimensional flow field. Geometry is input to the program in the form of equa- 
tion coefficients. ODPART is equipped for conic and polynominal equations (see 
B2.1, Detailed Guide for Input Data). During a one-dimensional calculation the 
local radius from the nozzle centerline to the nozzle wall is calculated as a func- 
tion of axial coordinate using the input equation coefficients. Program logic aro- 
matically switches from one equation to another as the maximum axial coordinate 
is reached for which a given equation is applicable. A maximum of 10 sets of 
equation coefficients may be input to define the nozzle boundary. 

Data must be input to define individual particle sizes, particle mass den- 
sities and the percentage of the total particle mass comprised by each particle. 
The total particle mass relative to the gas must also be input. These proper- 
ties define the physical characteristics of the particle phase. A maximum of 
six different particle species may be input. To define the particle thermody- 
namics, the user may input either a temperature- enthalpy table or thermody- 
namic properties for each particle species. A temperature-enthalpy table con- 
sists of tabulated values of particle temperature and the particle enthalpy at each 
temperature. Tne temperature range should include the temperature at which 
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the particle species changes phase from solid to liquid. A maximum of 10 tem- 
perature-enthalpy points may be input for each particle species. If the user 
chooses to input thermodynamic and phase change properties, particle temper- 
ature and enthalpy will be calculated using ideal gas relationships. The ther- 
modynamic and phase data required for inpu* are melting temperature, the en- 
thalpy of the solid and liquid phases at the melting temperature and the speci- 
fic heat at constant pressure for the solid and liquid phases. 

To define the drag charcteristics of the particle species, a table of Rey- 
nolds number and corresponding drag coefficients must be input. A maximum 
of 50 Reynolds number- drag coefficients points may be input. The user may in - 
put a particle drag tabl ' or use a drag table contained in the program. If the 
program control flag, IDRAG, is set to 0, the program will use the particle drag 
table cc itained internally (Kliegel drag data, see Ref. 7). If IDRAG is set to 1, 
the user must input a drag table. 

The numerous calculational and input options of the ODPART computer 
code are outlined in greater detail in Section B-2. 

B.2 DESCRIPTION OF PROGRAM INPUT 

This subsection contains a detailed description of the program input as 
follows: 


• Detailed input guide 

• Detailed description of the input FORTRAN 
symbols 
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Fig. B- 1 - Schematic of One-Dimensional Gas-Particle Solution 
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B. 2. 1 Detailed Guide for Input Data 


Card 1 - Program control flags Format: 1615 


Column 

Parameter 


Description 

1-5 

ICON(l) 

0 

Particle limiting streamlines will 
be calculated 



1 

Particles will be distributed on 
start line for a nozzle flowing full 

6-10 

ICON(2) 

0 

Straight start line (constant axial 
coordinate) 



1 

Normal start line (points calcula- 
ted at intersection of source stream 
lines with curved line at start line 
position) 

11-15 

ICON (3) 

Number of start line points 

16-20 

ICON(4) 

Number of upper boundary equations 
(maximum of 10) 

21-25 

ICON(9) 

0 

Use English units for nozzle geo- 
metry 

26-30 

ICON(IO) 

0 

Geometry input in inches 



1 

Geometry input in feet 

31-35 

ICON(ll) 

0 

Punch start line points 



1 

Do not punch start line points 

36-40 

ICON(12) 

0 

Points distributed on start line with 
sine function 



1 

Points distributed on start line at 
even radial increments 

41-43 

IC-ON(13) 

0 

Start line calculated for input to 
RAMP equilibrium code 



1 

Start line calculated for input to 
RAMP finite rate code 

Card Z - Calculational control parameters Format: 1615 

1-5 

ISPECS 

Number of particle species (maximum 
of 6) 

* 

6-10 

IDRAG 

0 

Use drag table coded in program 


N N imber of points to be input for new 
drag table 

- 

Kliegel drag data (Ref, 7) 
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Card 2 - (Con' t) 


Column 

Parameter 

Description 

11-15 

IDEL 

Interval at which properties of data 
surfaces are printed 

16-20 

ISAV 

Interval at which properties of data 
surfaces are stored for use in par- 
ticle limiting streamline calculations 


Card 3 - Initial conditions and problem limits Format: 7E10.6 


1-10 

XM1 

Axial coordinate of throat 

11-20 

XI 

Axial coordinate of initial data surface 

21-30 

DX 

Axial integration step 

31-40 

DXJ 

Throat radius 

H- 50 

UG 

Gas velocity on initial data surface. 
IF XSTART >0, UG = 0 

51-60 

UP 

Particle velocity on initial data surface. 
IF XSTART >0, UP = 0 

61-70 

XSTART 

Axial location of initial data surface if 
internal calculation of initial gas and 
particle velocities is desired 

Tl-80 

EM OUT 

Mach number at which calculations are 
to stop and start line to be punched 


Card 4 - Gas -particle equilibrium thermodynamic Format: 15, 5X, El 0.6, A3 
table (only one table required) 


1-5 

IVEQ 

Number of Mach numbers in the gas- 
particle equilibrium thermo table 

11-20 

HC 

Chamber enthalpy 

21-30 

UNITS 

ENG Input gas data with English units 
MKS Input gas data with Metric units 


Card 4 - Gas-particle equilibrium thermodynamic Format: BE10. 6 
data (repeat Card 4 IVEQ times in order 
of increasing Mach number) 

1*10 EQUIL(J, 1) Gas-particle equilibrium Mach number 

11-20 EQUIL(J, 2) Gas-particle equilibrium molecular 

weight 
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Card 4 - (Con' t) 


Column 

Parameter 

Description 

21-30 

EQUILtf, 3) 

Gas -particle equilibrium ratio of 
specific heats 

31-40 

EQUILtf, 4) 

Gas-particle equilibrium temperature 

41-50 

EQULL(J, 5) 

Gas -particle equilibrium pressure 

51-60 

EQUILtf, 6) 

Gas-particle equilibrium Prandtl 
number 

61-70 

EQUILtf, 7) 

Gas-particle equilibrium viscosity 

71-80 

EQUILtf, 8) 

Gas-particle equilibrium specific 
heat at constant pressure 

Card 5 - Control flags for gas phase thermodynamic Format: 215 
tables 

1-5 

me 

Number of enthalpy tables (maud- 
mum of 10) 

6-10 

LSA 

Number of entropy cuts for each en- 
thalpy table (maximum of 2) 

Card 6 - Total enthalpv of table 

Format: El 0.6 

1-10 

DELHC(M) 

Enthalpy from which gas was ex- 
panded 

Card 7 - Entropy of table and no. 
at this entropy 

Mach numbers 

1-10 

STAB(M, I) 

Entropy at which gas was expanded 

19-20 

IVTAB(M, I) 

Number of Mach numbers in thermo- 
dynamic table 

Card 8 - Gas phase thermodynamic data 

1-10 

GAStf, 1) 

Gas phase Mach number 

11-20 

GAStf, 2) 

Gas phase molecular weight 

21-30 

GAStf, 3) 

Gas phase ratio of specifi: heats 
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Card tf - (Con 1 1) 


Column 

Parameter 

Description 

31-40 

DUM 

Not currently used 

41-30 

DUM 

Not currently used 

51-60 

GAS(J, 4) 

Gas phase Prandti number 

61-70 

GAS(J, 5) 

Gas phase viscosity 

71-80 

GAS(J, 6) 

Gas phase C 


Note: Repeat Cards 6, 7 and 8 IHC times in order of increasing total 
enthalpy, DELHC(M) 


Card 9 - Coefficients of equations defining Format: II, 3X, II, 5X, 

nozzle geometry 6E10.6 

1 IWALL(K, 2) 1 Conic Equation 7 x / ? 

R = A[(B + CX + Da ) A/ + E] 

2 Polynomial Equation 

R = AX 4 + BX 3 + CX 2 + DX + E 


5 

ITRAN (K, 2) 

Not presently used 

11-20 

WALLCO 
(K, 1, 2) 

Coefficient A 

21-30 

WALLCO 
(K, 2, 2) 

Coefficient B 

31-40 

WALLCO 
(K, 3, 2) 

Coefficient C 

41-50 

WALLCO 
(K, 4, 2) 

Coefficient D 

51-60 

WALLCO 
(K, 5, 2) 

Coefficient E 

61-70 

WALLCO 
(K, 6, 2) 

Maximum value of X applicable to 
equation 


Note: Repeat Card 9 sufficient number of times to define the nozzle 

inlet, throat and divergent exit 


B - 10 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



Card 10 


Column Parameter Description 

1-10 XMASSP Ratio of particle total mass flowrate to 

gas mass flowrate 


Card 1 1 - Percentage of particle mass flow Format: 6E10.6 

contributed by each particle 

I- 10 PERTG(l) Ratio of particle 1 mass flowrate to 

total particle mass flowrate 

II- 20 PERTG(2) Ratio of particle 2 mass flowrate to 

• * total particle mass flowrate 

• • • 

• • • 

PERTG • 

(SPECS) 


Card 12 — Radius of each particle Format: 6E10.6 

1-10 PSP(2, 1) Radius of particle 1, microns 


PSP(2, SPECS) Radius of particle SPECS, microns 


Card 13 — Mass density of each particle Format: 6E10.6 

1-10 PSP(1, 1) Mass density of particle 1 



PSP(1, SPECS) 


Card 14 — Particle drag data (Use only if IDRAG>0) Format: 8E10.6 

I- 10 AREY(l) Reynolds number used as independent 

variable in drag table 

II- 20 BREY(l) Drag coefficient corresponding to 

above Reynolds number (Cq/Cjj stokes^ 
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Card 14 — (Con 1 1) 


Column 

Parameter 

Description 

21-30 

AREY (2) 


31-40 

BREY(2) 


• 

• 


• 

• 


• 

• 



Use as many cards as necessary to input all necessary drag data (50 max 
entries) 


Card 1 5 — Particle data control variables Format: 4A6, B, A6 

I -24 ALPHA Particle name 

25-27 NOCUTS Not used 

23-33 UNIT ENG Data input in English units 

MKS Data input in Metric units 

(left adjusted) 


Card 16 — Thermodynamic table to be used by each Format: 1012 

particle species 

1-2 JTEM(l) Temperature- enthalpy table to be used 

• • for particle 1 


JTEM(ISPECS) Temperature-enthalpy table to be used 
for particle ISPECS 

Card 17 — Number data points in particle thermo- Format: B 

dynamic table 

1-3 NPTM(I) Number of temperature-enthalpy data 

points for this particle. If equal to 1, 
input liquid and solid heat capacities. 
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Card 18 — Particle thermodynamic data 


Format: 7E10.6 


Column 

Parameter 

Description 

1-10 

TM(I) 

Melting point temperature of particle I 

11-20 

Ho(I) 

Enthalpy of solid phase of particle I at 
melting point temperature 

21-30 

HM(I) 

Enthalpy of liquid phase of particle I at 
melting point temperature 

If NPTfc (I) = 1 

, use the following 

format 

31-40 

APHO(l, 1, I) 

Heat capacity of liquid phase of particle I 
(Btu/lbm °R or cal/gm-°K) 

41-50 

APHO(l, 2, I) 

Heat capacity of solid phase of particle I 
(Btu/lbm"°R or cal/gm-K) 

If NPTM(I)>1, 

use following format 

31-40 

APHO(l, 1, I) 

Temperature for T-H table for particle I 
(°R or °K) 

41-50 

APHO(l, 2, I) 

Enthalpy for T-H table for particle I (Btu/ 
lbm or cal/gm) 

51-60 

APHO(2, 1, I) 

Second temperature in T-H table for par- 
tide I (°R or °K) 

61-70 

APHO(2, 2, I) 

Second enthalpy in T-H table for particle I 
(Btu/lbm or cal/ gm) 


The above format (APHO(J, 1, I), APHO(J, 2, I)) is continued on successive 
cards of format 7E10. 6 for J = 1, Z, . . . NPTM(I). 


B. Z. Z Input FORTRAN Symbols 


Symbol 

Description 

Units 

ICON(l) 

This flag designates type of particulate start 
line to be calculated and punched. If ICON(l) 
= 0, particle trajectories and their intersec- 
tion with the gas start line are calculated. If 
ICON(l) = 1, the particles are evenly distri- 
buted across the gas start line. 

n/a 

ICON(2) 

This flag designates type of gas start line to 
be calculated and punched. If ICON(2) = 0, a 
straight start line is calculated. If ICON(2) 

= 1, a normal start line is calculated. 

n/a 


B- 13 



LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



Symbol 

Description 

Units 

ICON(3) 

ICON(3) is the number of start line points 
to be generated. A maximum of 50 is per- 
mitted. 

n/a 

ICON(4) 

ICON(4) is the number of upper boundary 
equations defining the nozzle geometry. 

n/a 

ICON( 9) 

ICON(9) specifies the units to be used* 
If ICON(9) = 0, English units are used. 
If ICON(9) = 1, Mks units are used. 

n/a 

ICON(IO) 

This flag designates the dimensions of the 
input geometry and geometric control para- 
meters if English units are used. If ICON(IO) 
= 0, geometry is input in inches. If ICON(IO) 
= 1, geometry is input in feet. 

n/a 

ISPECS 

Number of particle species to be input. A 
maximum of 10 is permitted. 

n/a 

IDRAG 

IDRAG controls which particle drag data is 
used. If IDRAG = 0, the particle drag table 
coded in the program will be used. If a drag 

table is input, IDRAG is number of points to 
be read. 

n/a 

IDEL 

IDEL is the interval at which properties of 
data surfaces are printed. If IDEL = 1, the 
properties of each data surface are printed. 
If IDEL = 10, the properties of every tenth 
data surface are printed, etc. 

n/a 

ISAV 

ISAV is the interval at which properties of 
data surfaces are stored for use in particle 
trajectory calculations. A maximum of 300 
data surfaces may be stored. 

n/a 

XM1 

Axial coordinate of nozzle throat 

in or ft 

XI 

Axial coordinate of initial data surface 

in or ft 

DX 

Integration step between successive data 
surface calculations 

in or ft 

DXJ 

Radius of the nozzle throat 

in o r ft 

UG 

This variable is the gas velocity on the initial 
data surface. If the velocity is unknown, UG 
should be set equal to zero and the program 
will calculate an initial guess for UG. 

ft/ sec 

UP 

This variable is the particle velocity on the 
initial data surface. If UP is set equal to 
zero, the program will set UP = UG and as- 
sume the gas and particles are in equilibrium 
on the initial data surface. 

ft/ sec 


B - 14 


LOCKHEEO • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



Symbol 

Description 

Unite 

XSTART 

XSTART is the axial location of the initial 
data surface if internal calculation of ini- 
tial gas and particle velocities is desired. 

in or ft 

EMOUT 

This is the Mach number at which the inte- 
gration of the gas-particle equations is to 
stop and a start line calculated and punched. 

n/a 

IV EQ 

This flag designates the number of Mach num- 
ber entries in the gas-particle equilibrium ther- 
modynamic table. A maximum of 13 is permit- 
ted. 

n/a 

HC 

This is the chamber enthalpy in units consis- 
tent with the thermodynamic tables. 

ft 2 / sec c 

-> or 2 
M 2 /sec 

EQUIL.(J, 1) 

Mach number entry in the gas -particle equi- 
librium thermodynamic table. 

n/a 

EQUIL(J, 2) 

Molecular weight entry in the gas -particle 
equilibrium thermodynamic table 

lbm./lbm. mole 
or 

gm/grr*. mole 

EQUU.(J, 3) 

Ratio of specific heats entry in the gas- 
particle equilibrium thermodynamic table 

n/a 

EQUIL(J, 4) 

Temperature entry in the gas-particle equi- 
librium thermodynamic table 

°R 

or 

°K 

EQULL(J, 5) 

Pressure entry in the gas-particle equilibrium 
thermodynamic table 

atm. 

EQUIL(J, 6) 

Prandtl number entry in the gas-particle equi- 
librium thermodynamic table 

n/a 

EQUIL(J, 7) 

Viscosity entry in the gas-particle equilibrium 
thermodynamic table 

lbf. sec 

ft 2 

or 

Kg- _m 

sec 

EQUIL(J, 8) 

Gas specific heat at constant pressure entry 
in the gas-particle equilibrium thermodynam- 
ic table 

ft 2 

2 o D 
sec . R 

IHC 

Number of enthalpy tables in the gas phase 
thermodynamic data A maximum of 10 is 
permitted. 

n/a 

ISA 

Number of entropy cuts for each enthalpy 
table. A maximum of 2 is permitted. 

n/a 
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Symbol 

Description 

Units 

DELHC(M) 

This variable is the total enthalpy from 
which the gas was expanded. 

Btu/lbm 

or 

cal/ gm 

STAB(M, I) 

This variable is the entropy at which the 
gas was expanded 

Btu/lbm - °R 
or 

cal/ gm - K 

IVTAB(M, I) 

This flag designates the number of Mach 
number entries in the gas phase thermodyna- 
mic table. A maximum of 1 3 is permitted. 

n/a 

GAS(J, 1) 

Mach number entry in the gas phase thermo- 
dynamic table. 

n/a 

GAS{J, 2) 

Molecular weight entry in the gas phase 
thermodynamic table 

lbm/lbm. mole 
or 

gm/ gm. mole 

GAS(J, 3) 

Ratio of specific heats in the gas phase 
thermodynamic table. 

n/a 

GAS ( J , 7) 

Temperature entry in the gas phase thermo- 
dynamic table. 

°R 

or 

o,_ 

GAS(J, 4) 

Prandtl number entry in the gas phase thermo- 
dynamic table. 

K 

n/a 

GAS(J, 5) 

Viscosity entry in the gas phase thermody- 
namic table. 

lbf. sec 

— z — 

ft 

or 

kg. M 
sec 

GAS(J, 6) 

Gas specific heat at constant pressure entry 
in the gas phase thermodynamic table. For 
ideal gas thermodynamics it is the non-di- 
mensional temperature exponent. 

ft 2 ft 2 

2 o 
sec - R 

or 

M 2 

sec 2 - °K 

IWAL(K, 2) 

This flag designates the type of equation de- 
fining the nozzle geometry. If IWALL(K, 2) 
= 1, a conic equation of the form, R = A[(B + 

CX + DX 2 ) 1 ^ 2 + E] is input. If IWALL1K, 2) 
= 2, a polynomial equation of the form, R = 

AX 4 + BX 3 + CX 2 + DX + E is input. 

n/a 
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Symbol 
ITRAN(K, 2) 
WALLCO(K, 1,2) 

WALLCO(K, 2, 2) 

WALLCO(K, 3, 2) 

WALLCO(K,4, 2) 

WALLCO(K, 5, 2) 

WALLCO(K, 6, 2) 

XMASSP 

PERTG(J) 

PSP(2, J) 

PSP(1, J) 

ALPHA 

NOCUTS 

UNIT 

JTEM(J) 

NPTM(J) 


Description 
Not presently used. 

This variable corresponds to the coef£icient 
A in the conic or polynomial nozzle equations. 


This variable corresponds to the coefficient 
B in the conic or polynomial nozzle equations. 

This variable corresponds to the coefficient 
C in the conic or polynomial nozzle equations. 

This variable corresponds to the coefficient 
D in the conic or polynomial nozzle equations. 

This variable corresponds to the coefficient 
E in the conic or polynomial nozzle equations. 

Maximum axial coordinate for which a parti- 
cular nozzle equation is applicable. 

XMASSP is the ratio of total particle mass 
flow rate to gas mass flow rate 

This is the ratio of the Jth particle mass 
flow rate to total particle mass flow rate 

PSP(2, J) is the radius of the Jth particle 

PSP(1, J) is the mass density of the Jth 
particle 


ALPHA is the particle specie name is used 
for printout purposes. 

Not presently used 

UNIT specifies the units in which the particle 
temperature and enthalpy data will be input. 

If UNIT = ENG, data will be input in English 
units. If UNIT = MKS, 'data will be input in 
MKS units. 

JTEM(J) specifies the temperature- enthalpy 
table to be used for the Jth particle 

Number of entries in temperature-enthalpv 
table for a particular particle. If NPTM(J) 

= 1, input the liquid and solid heat capacities. 
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Units 

N/A 

Equation coef- 
ficients in the 
WALLCO array 
should be devel- 
oped in units con- 
sistent with ICON 
ICON(IO) 


n/a 

n/a 

Microns 

lbm/ ft^ 
or 

gm/ cm^ 

n/a 

n/a 

n/a 


n/a 

n/a 
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Symbol 

inscription 

Units 

TM(J) 

Melting point temperature of Jth particle 

°R 

or 

HS(J) 

Enthalpy of solid phase of Jth particle at 
the melting point temperature 

°K 

Btu/lbm 

or 

cal/ gm 

HM( J) 

Enthalpy of liquid phase of Jth p i tide at 
the melting point temperature 

Btu/lbm 

or 

cal/ gm 

APHO(l, 1, J) 

Heat capacity of the liquid phase of the Jth 
particle if NPTM(J) = 1. First temperature 
entry in the Jth particle temperature- enthalpy 
table if NPTM(J)> 1. 

sec detailed 
input guide 

APHO( 1, 2, J) 

Heat capacity of the solid phase of the Jth 
particle if NPTM(J) = 1. First enthalpy entry 
in the Jth particle temperature-enthalpy table 
if NPTM(J) > 1. 

see detailed 
input guide 

APHO(z, 1, J) 

Second temperature entry in the Jth particle 
temperature-enthalpy table if NPTM(J)> 1. 

°R 
or 
o „ 

APH(z, z, J) 

Second enthalpy entry in the Jth particle 
temperature-enthalpy table if NPTM(J)>1. 

K 

Btu/lbm 

or 

cal/ gm 


B. 3 DESCRIPTION OF PROGRAM OUTPUT 

This subsection contains a description of the output scheme utilized by 
the program and a description of the error messages printed out by the pro- 
gram. 

B. 3. 1 Description of Program Data Output 

The following guide consist^ numbered flags on sample printout sheets 
(pages B-25 through B-32) which respond to numbered comments listed below, 

1 Program title 

2 Control flags input by user 

3 THROAT X COORDINATE Axial coordinate of r. ~zle throat, 

ft or in 
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4 

AXIAL COORDINATE AT 
THE THROAT ENTRANCE 

Axial coordinate of nozzle inlet, 
ft or in 

5 

DELTA X 

Axial integration step, ft or in 


♦ 

R 

Nozzle throat radius, ft or in 

7 

MIXTURE VELOCITY AT 

Velocity at nozzle inlet input by user, 
ft/sec 

8 

CHAMBER El JHALPY 

Chamber enthalpy for wni^u gas-particle 

thermodynamics were calculated, 

f 2 / 2 x .2 / 2 
ft /sec or M /sec 

9 

VELOCITY 

Velocity for which gas-particle thermo- 
dynamics were calculated, ft/sec 

10 

R 

Gas constant, ft^/sec^ - °R or M^/sec^- 

11 

GAMMA 

Ratio of specific heats 

12 

TEMPERATURE 

Static temperature, °R or °K 

13 

PRESSURE 

Static pressure, psfa 

14 

PRANDTL NUMBER 

Prandtl number 

15 

VISCOSITY 

Viscosity, poise 

lb 

CP 

Specific heat at constant pressure, 
ft^/sec^ - °R or M^/sec^ - °K 

17 

ENTHALPY 

Chamber enthalpy fcr which gas phase 

thermodynamics were calculated, 

.2/ 2 xx 2 2 
ft / sec or M / ^ec 

18 

ENTROPY 

Chamber entropy for which gas phase 
thermodynamics were calculated, ft*V 
sec - R or M /sec - K 

19 

VELOCITY 

Velocity for which gas-particle thermo- 
dynamics were calculated. 

20 

R 

Ga onstant, ft^/sec^ - °R or M^/sec^ - 

21 

GAMMA 

Ratio of specific heats 

22 

PRANDTL NUMBER 

Prandtl number 

23 

VISCOSITY 

Viscosity, poise 

24 

EXPONENT 

Viscosity exponent 

25 

TYPE 

* ype of nozzle boundary equation 
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26 

ITRANS 

Gas dynamic condition existing at the 
end of the current nozzle boundary 
equations (not currently used) 

27 

A 

Nozzle boundary equation coefficient 

28 

B 

Nozzle boundary equation coefficient 

29 

C 

Nozzle boundary equation coefficient 

30 

D 

Nozzle boundary equation coefficient 

31 

E 

Nozzle boundary equation coefficient 

32 

MAX 

Maximum axial coordinate for which 
nozzle boundary equation is appicable 

33 

TOTAL PARTICLE WEIGHT 
RELATIVE TO THE GAS IS 

Ratio of total particle mass flow rate 
to gas mass flow rate 

34 

I 

Particle number 

35 

RADIUS 

Particle radius, microns 

36 

PERCENT 

Ratio of Jth particle mass flow rate to 
total particle mass flow rate 

37 

DENSITY 

Particle mass density, lbm/ft' 

38 

VELOCITY 

Particle velocity at first axial station, 
ft/ sec 

39 

TMELT 

Melting temperature of particle species, 
°R or °K 

40 

HSOLID 

Enthalpy of solid particle at melting tern- 
perature, ft /sec or M /sec 

41 

HLIQUID 

Enthalpy of liquid particle at melting tern 

r Z / 2 7/ 2 

perature, ft /sec or M /sec 

42 

CPMELT 

Heat capacity of liquid particle at melt- 
ing temperature, ft^/sec^ - °R or M^/ 

2 

sec - K 

43 

CPSOLID 

Heat capacity of solid particle at melt- 

2 2 o 2 

ing temperature, ft /sec - R or M / 
2 

sec - K 

44 

I 

Point number 

45 

RE 

Reynolds number 

46 

DRAG COEF 

Drag Coefficient (Cp/Cp ) 

^Stokes 
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4? 

AXIAL STATION X 

Axial coordinate of data line, ft or in 

48 

AREA RATIO A/A * 

Area ratio at the axial location of the 
data line 

49 

VELOCITY 

Gas velocity, ft/ sec 

50 

R 

2 2 q 

Gas constant, ft /sec - R 

51 

ENTHALPY 

2 2 

Gas static enthalpy, ft / sec 

52 

GAMMA 

Ratio of specific heats 

53 

TEMPERATURE 

Gas static terrperature, °R 

54 

PRANDTL NUMBER 

Prandtl number 

55 

MACH NUMBER 

Mach number 

56 

VISCOSITY 

Viscosity, poise 

57 

PRESSURE 

Gas static pressure, psia 

58 

CP 

Specific heat at constant pressure, 
ft^/ sec^ - °R 

59 

DENSITY 

Gas density, slugs/ ft 

60 

TOTAL ENTHALPY 

2 2 

Gas total enthalpy, ft /sec 

61 

ENTROPY 

Entropy, ft^/sec^ - °R 

62 

PARTICLE 

Particle number 

63 

VELOCITY 

Particle velocity, ft/ sec 

64 

ENTHALPY 

2 / 2 

Particle enthalpy, ft /sec 

65 

TEMPERATURE 

Particle temperature, °R 

66 

MIXTURE MASS FLOW 

Percent change in mass flow from the 
initial data line to the calculated start 
line 

67 

GAS MOMENTUM 

Percent change in gas momentum from 
the initial data line to the calculated 
start li^e 

58 

PARTICLE MOMENTUM 

Percent change in particle momentum 
from the initial data line to the calcu- 
lated start line 

69 

MIXTURE MOMENTUM 

Percent change in mixture momentum 
from the initial data line to the calcu- 


lated start line 
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70 

TOTAL ENTHALPY 

Fercent change in total enthalpy from 
the initial data line to the calculated 
start line 

71 

FXG 

Integral of force parallel to the nozzle 
centerline due to the gas, lbf 

72 

FXP 

Integral of force parallel to the nozzle 
centerline due to the particles, lbf 

73 

FXM 

Total integral of force parallel to the 
nozzle centerline, lbf 

74 

ISP 

One-dimensional ISP, sec 

75 


Particle number 

76 

RADIAL COORDINATE 

Radial distance from nozzle centerline, ft 

77 

AXIAL COORDINATE 

Axial lot ation, ft 

78 

AXIAL VELOCITY 

Particle axial velocity, ft/ sec 

79 

RADIAL VELOCITY 

Particle radial velocity, ft/sec 

80 

ANGLE 

Angle from the horizontal of the par- 
ticle trajectory, degrees 

81 

ENTHALPY 

2 2 

Particle enthalpy, ft / sec 

82 

NOZZLE RADIUS 

Nozzle radius, ft 

83 


Particle number 

84 

RADIAL COORDINATE 

Radial coordinate of particle limiting 
streamline intersection with calculated 
start line, ft 

85 

AXIAL COORDINATE 

Axial coordinate of particle limiting 
streamline intersection with calculated 
start line, ft 

86 

AXIAL VELOCITY 

Particle axial velocity at intersection 
of particle limiting streamline with cal- 
culated start line, ft/ sec 

87 

RADIAL VELOCITY 

Particle radial velocity at intersection of 
particle limiting streamline with calculated 
start line, ft/ sec 

88 

ANGLE 

Angle from the horizontal of the particle 
trajectory at intersection of particle 
limiting streamline with calculated start 
line, degrees 

89 

ENTHALPY 

Particle enthalpy at intersection of par- 
ticle limiting streamline with calculated 
2 / 2 

start line, ft / sec 
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90 

NOZZLE RADIUS 

Nozzle radius, ft 

91 

R 

Radial coordinate of start line point, ft 

92 

X 

Axial coordinate of start line point, ft 

93 

MACH NUMBER 

Mach number at start line point 

94 

THETA 

Gas flew angle at start line point, degrees 

95 

ENTROPY 

2 2 o 

Entropy at start line point, ft / sec - R 

96 

MACH ANGLE 

Mach angle at start line point, degrees 

97 

TOTAL ENTHALPY 

Total enthalpy at start line point, 

,2/ 2 
ft /sec 

98 

PARTICLE 

Particle number 

99 

I 

Start line point number 

100 

AXIAL VELOCITY 

Particle axial velocity at start line point, 
ft/ sec 

101 

RADIAL VELOCITY 

Particle radial velocity at start line 
point, ft/ sec 

102 

THETA 

Angle from the horizontal of the parti- 
cle trajectory at the start line point, 
degrees 

103 

ENTHALPY 

Particle enthalpy at the start line point, 

.Z i 2 

ft / sec 

104 

DENSITY 

Particle denisty at the start line point, 
Ibm/ ft^ 
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Note: These cards were input to a Univac 1108 multiprocessor system to generate the 

sample output which follows. 
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RUN CONTROL RuRRHCURS 
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particle trajectory calculations 
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B.3.2 Description of Program Srror Message Output 

1. THE MIXTURE START CONDITION WAS NOT FOUND 
IN 300 STEPS 

This statement is printed if the program did not reach 
the desired startline Mach number after 300 axial sta- 
tioi s were stored. 

2. THE GAS PARTICLE FLOW SOLUTION AT X= AND A= 

WILL NOT CONVERGE 

This statement is printed if the program was unable to 
converge on a solution after 100 iterations at a parti- 
cular axial location. The run is terminated following 
this message. 

3. PARTICLE SOLUTION WILL NOT CONVERGE AT X= 

This message is printed if the particle trajectory tracing 
routine is unable to converge on a solution after 50 iter- 
ations at a particular axial station. The run is continued 
and a startline generated for the last axial station at 
which the solution converged. 

4. ITSUB WNC 

This rm ssage originates in subroutine ITSUB which pro- 
vides iteration control for a given function. The message 
indicates that a solution was not obtained after 100 itera- 
tions. 

5. CHAM3ER ENTHALPY NOT FOUND IN SUBIN 

This message originates in subroutine SUBIN and indicates 
that the routine was unable to find a gas phase thermodynamic 
table expanded irom the input chamber enthalpy. 
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